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We study the dynamical generation of masses for fundamental fermions in quenched quantum 
electrodynamics in the presence of weak magnetic fields using Schwinger-Dyson equations. Contrary 
to the case where the magnetic field is strong, in the weak field limit the coupling should exceed 
certain critical value in order for the generation of masses to take place, just as in the case where no 
magnetic field is present. The weak field limit is defined as eB -C m(0) 2 , where m(0) is the value 
of the dynamically generated mass in the absence of the field. We carry out a numerical analysis to 
study the magnetic field dependence of the mass function above critical coupling and show that in 
this regime the dynamically generated mass and the chiral condensate for the lowest Landau level 
increase proportionally to (eB) 2 . 
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It is well known that in QED, fermions can acquire 
masses through self interactions without the need of a 
nonzero bare mass. This phenomenon, known as dy- 
namical mass generation (DMG), happens above a cer- 
tain critical value of the coupling and its description can 
only be carried out in terms of non-perturbative treat- 
ments. Schwinger-Dyson Equations (SDEs) provide a 
natural platform to study DMG. In the quenched version 
of QED, a favorite starting point is to make an ansatz 
for the fermion-photon vertex and then study the fermion 
propagator equation in its decoupled form. It is also well 
known that in the presence of strong magnetic fields, it is 
possible to generate fermion masses for any value of the 
coupling. This phenomenon has been given the name of 
magnetic catalysis 0- [3- 0- 01 ■ Non-perturbative aspects 
of dynamical mass generation in the presence of weak 
magnetic fields have earlier been considered in the con- 
text of Nambu-Jona-Lasinio (NJL) model QCD Q 
and (2+I)-dimensional QED jj- In the context of QED4, 
the only work to our knowledge, is that of Kikuchi and 
Ng 8]. However, they concentrate mainly on the be- 
haviour of the critical coupling in the presence of the 
weak magnetic fields. In this paper, we undertake the 
study of the weak field dependence of the dynamically 
generated mass and the chiral condensate in QED in the 
rainbow truncation of SDE. 

The SDE for the fermion propagator without external 
fields (in vacuum) in the rainbow approximation is 

s-\ P ) = s^-\p) -^Jd A k -fs F (k)-fA$(q),W 

where q — k — p and a = e 2 /(4ir) is the electromag- 
netic coupling constant. In this expression A^(q) is 
the bare photon propagator, which in covariant gauges 
is written as A$(q) = - (g^ + (£ - l)q^qv/q 2 ) /q 2 , £ 
being the usual covariant gauge parameter. We write the 
full fermion propagator as Sf(p) = F(p 2 )/(p'— M(p 2 )). 



F(p 2 ) is referred to as the wave function renormalization 
and M(p 2 ) as the mass function. In the Landau gauge 
(£ = 0) F(p 2 ) — I and the mass function has nontriv- 
ial solutions for values of the coupling above the critical 
value a c — ir/3. In the presence of external fields, SDEs 
have been a subject of study already for some time, see 
for example Ref. 0- 

When the magnetic field is strong, Landau levels are 
separated from each other by an amount ~ v 'eB in such 
a way that for any value of the coupling a, only the lowest 
Landau level (LLL) contributes to the DMG HH 0] . 
However, in the case of weak external magnetic fields, 
Landau levels are close to each other and hence all contri- 
butions should be taken into account, which adds consid- 
erably to the complexity of the problem, as emphasized 
also in the fourth article of reference [fj . 

The presence of the field breaks Lorentz invariance. 
Consequently, a simple Fourier transform on a single mo- 
mentum variable is not possible. Nevertheless, it has 
been shown 0] that the mass operator in the presence 
of an electromagnetic field can be written as a combina- 
tion of the structures 

7^ , a^F, u , (F^IT) 2 , 75 (2) 

which commute with the operator (7 • n) 2 , where 
LL, = idp-eA?, = d„A°* -d u A°f, F^ = 

i e fifXrp^^ o ^ v = 7„]/2 and A° xt is the external vec- 
tor potential. We take A cxt = ff(0, -y/2, x/2 0) which 
describes a constant magnetic field B = Bz [TI| . 

In order to find a diagonal representation for the mass 
operator, we thus need to find the eigenfunctions ip pa of 
the operator (7 • II) 2 , namely 

(7 • U) 2 1p pa U ax = p 2 1pp a Ucr X , (3) 

where taken as the eigenspinors of S 3 and 75. 

We work in cylindrical coordinates r = (r, <f>, z) and in 
the chiral representation of the 7— matrices where S3 and 



2 



75 are both diagonal with eigenvalues a = ±1 and x = 
±1. The normalized cigenfunctions ip p(T (see for example 
Ref. 0|) are gi ven by 



vM (4) 



to find the explicit form for the function Q in momentum 
space, which is given by 



Q(k) 



1 



7- k- E(fc) 



(9) 



where N = ^27/(2tt) 3 , p = 7r 2 , 7 = eS/2, p 2 
p 2 — 2eBn and 



I-u(p) 



-p/2 p {n-s)/2 L n- s{j ^ 



(5) 



are the Laguerre functions |13| with the quantum num- 
bers n,l,s related by n — l+s. Since the problem involves 
only a magnetic field, the solutions do not depend on the 
eigenvalues \- 

The solutions can be conveniently arranged in a matrix 
form 



(6) 



cr=±l 



where A(cr) = diag <5<r-i> 6<ri, #<r-i} is a 4 x 4 matrix 
and x = (t, r). 

The matrix in Eq. @ is used to rotate the two-point 
fermion Green's function between coordinate, G{x, y) 
and momentum spaces, Q(k), as 



G(x,y) = J d 2 k^ k (x)g(k)^ k (y), 



(7) 



where k\\ = (fc o ,0,0, fc 3 ) and * fe = 7°*j c 7°. The above 
expression can be substituted into the equation relating 
the two-point fermion Green's function and the mass op- 
erator M(x, y) in coordinate space, namely 



In order to arrive at this equation, we have used the 
completeness of the functions ipka- expressed in terms of 
*fc as 



£ I d z k^ k {x)^ k {y) = 8\x-y) (1.0! 

"fc.Sfc ' 



along with the properties 



7-n(x)#fc(a;) = # fc (x)(7-fc) 

d A x'M{x,x')^ k {x') = * k (x)E(k), (11) 

and the definition of the mass operator S(fc) in momen- 
tum space 



M(k,k') = J d^xdtx'VkWMfax'^k'ix') 

= 5 nknk ,5 SkSk ,5 2 m-k[)Y.{k). (12) 

With the aid of Eqs. (|51)- (I12H it is now straightforward 
to transform the SDE for the mass operator in the rain- 
bow approximation from coordinate space, namely, 

M[x, x') = -ie 2 ^G(x, x')i»D$ (x - x') , (13) 



7 • H(x)G(x, y) — Id x'M(x, x')G(x' , y) = S i (x — y) (8) to momentum space, which now reads 



5n p n p ,5 SpSp ,5' 2 {p\\ -pj|)E(p) = -ie 2 I d 4 xd 4 x' 



d 2 k^ p {xh^ k (x) 



D^(x-x') 
7-fc-E(fc) 



Vkix'W^ix'), (14) 



where the bare photon propagator is 

/jA —iq-(x-x') 

Having considered the dependence of the mass func- 
tion S(fc) on the structures in Eq. J2J, its remain- 
ing, most general form can be written as fi — E(fc) = 
T^ 1 {k) [ft — M.(k)] , where, as in the case of vacuum, 
F(k) and A4(k) are called the wave function renormal- 
ization and mass functions, respectively, in the presence 



of the field. We work in the Landau gauge (£ = 0) where 
we know that for vacuum F = 1. Since we aim at a de- 
scription for small magnetic field strengths, we naturally 
expect T ~ 1 in the Landau gauge. Furthermore, let us 
work with the ansatz that M(k) is proportional to the 
unit matrix. Weakness of the magnetic field also implies 
that the bare vertex is a reasonable choice in the sense 
that Ward Identity is satisfied in the Landau gauge up 
to a correction connected with the mass function A4(k). 
This correction might be expected to be small because 
M(k) ~ M{p) for small values of momenta as the mass 
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function is practically a constant, and it falls off sharply 
as 1/k for large momenta. 

The self-consistent equation for the mass function is 
obtained by considering the diagonal part (n p — n p i s p — 
s p > ) and taking the trace of Eq. lfHI) . The integrals over 



x and x 1 in Eq. (|14|l are readily performed. Having set 
rip = n p i s p = s p >, the integral over x is just the complex 
conjugate of the one over x' . The first one is found from 
the expression 



47 J 



(l k -l P )/2-(a k -a p )/4 



e ^/^^-(^Va ( ^ /47)iS -.» ( ^ /47)j(16) 



where = (0,91,52,0). Equation (|16() is real and thus, 
integrating over x and x' in Eq. I|14f> . results in the square 
of the right-hand side of Eq. (|16[) . The presence of the 



delta function in this last equation, allows easy integra- 
tion over fen . Gathering the above described elements, 
we obtain self-consistent equation for the mass function 



M(p\\,n p ) 



E E 



d 4 q e -9i/ 2 7 



M((p- q)\\,n k ) 



2 \ lk—lp- 



4^ 



i" fc_ ; p ; + T7^(9i/47) 



Wr Sfc (gi/47)] 2 , (17) 



where in the notation for the mass function we have em- 
phasized the breakdown of Lorentz invariance. We ex- 
pect that M{{p — q)\\,rik) should be independent of Sfe 
since the energy only depends on the principal quantum 
number n k . Furthermore we assume that M,((p— q)\\,nk) 
is a slowly varying function of rife and thus make the ap- 
proximation M.((p — q)\\, rife) ~ M((p — q)\\,rik — 0). For 
consistency we consider the case n p = 0. Hereafter, we 
employ the more convenient notation M{ku,nk — 0) = 
M.{k\\) for generic arguments of the mass function. With 
these considerations the sum over s k can be computed by 



means of the result in Ref. 14] . It is worth mentioning 
that after summing over Sk, the resulting equation is the 
same as Eq. (50) in Ref. [2| when considering the case 
rife = 0, which corresponds to the strong field limit. 

In the situation where the magnetic field is weak, we 
expand [(p—q)^ — 2eBnk — M 2 ((p—q)\\)]~ 1 asageometric 
series in powers of eB. The remainin g su m over n k can 
be performed also by resorting to Ref. [lj] yielding, after 
a Wick rotation 



M(p» 



a 

47T 3 



-<?)||) 



9 2 [(p-?)jj 
6(6 - q\/q 2 )q\ 



<4 



M^ip-q^)] 



[(P-«)| 
36<zl 



■9i 



[(?■ 



q)l+ql+M 2 ((p-q) l \)] 3 [(p - «)fi + + 



A4 2 ((p-?)||)] 4 



-M 2 ((P-?)||)] 2 
(eS) 2 l , 



(18) 



keeping only the lowest order contribution in e73. Notice disappears on carrying out the sum over s/.. 

that, as expected, the s p dependence of the mass function . 

Solving the above equation numerically is still not 
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FIG. 1: Mass function in the LLL for different values of the 
weak external magnetic field for a — 6.5a c . 
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FIG. 2: Magnetic contribution to the dynamically generated 
mass (upper graph) and condensate (lower graph) in the LLL 
as a function of (eB) 2 for a — 6.5a c - 

trivial, owing to the fact that the unknown function 
■M((p — q)\\) within the integral is Lorentz non- invariant. 
However, we can always expand it out in powers of (eB) 2 . 
Therefore, M((p — q)\\) — Mo(p — q) + (eB) 2 M\, where 
M. i is responsible for breaking the Lorentz invariance of 
M (p-q). 

Consistently, we carry out the same expansion on the 
left hand side of Eq. H18|) . As the Lorentz invariance 
should be restored for the leading terms, we justifiably 
complete the momenta to achieve the same. This filters 
out the vacuum result. To calculate the magnetic field 
effect, we solve the integral equation for A^i obtained by 
comparing powers of (eB) 2 . The results for A4(p\\/A) in 
the LLL are depicted in Fig. ^ scaled by the ultraviolet 
cut-off A. Note that the results have been shown for 
a value of the coupling significantly above the critical 
value only for the magnetic field dependence to stand 
out. The same qualitative behaviour persists even in the 
immediate vicinity of the critical coupling, [l5| . 

The dot-dashed line corresponds to the vacuum. The 



effect of the external field is to increase the dynamically 
generated mass, preserving the qualitative features of the 
mass function profile. We have numerically verified that 
the critical value of the coupling is independent of the 
strength of the magnetic field, which is consistent with 
the findings of Ref . Q . To see the magnetic contribution 
to the dynamically generated mass, we show in Fig.[21the 
difference m(eB) — m(0), as a function of (eB) 2 , where 
m is the dynamical fermion mass, namely, m = A4(0). 
Notice that this difference grows linearly with (eB) 2 . 
This is the same behaviour as was observed for the NJL 
model in p|. We also evaluate the condensate defined 
as (ipil)) = i Tr G(x,x). In the weak field limit, spacing 
between Landau levels becomes small. To compute the 
condensate, the sum over these levels, can be carried out 
by replacing ^ n with the integral / d 2 k±/(2-K eB), along 
with the substitution 2eBn — > k\ . The weak field contri- 
bution to the condensate also turns out to be quadratic 
and with the above mentioned substitutions, it owes it- 
self entirely to the non Lorentz invariant piece of the mass 
function in our computational set up. Its behavior as a 
function of (eB) 2 is also shown in Fig. [21 

In summary, we have shown that for a > a c the dy- 
namically generated mass increases quadratically with 
the magnetic field strength. As compared to the strong 
field case, this is a four-fold dependence on the mag- 
netic field P, S S 01- Such a dependence is similar 
to that found in [5| for the NJL model. In the super- 
critical phase of QED, i.e, a > a c , the gap equation 
of Ref. y| can be solved in the linearized approxima- 
tion, giving the same quadratic dependence as we have 
demonstrated through explicit numerical evaluation of 
the SDE in our setup [17|- It is ineteresting to note 
that Farakos et. al. [7| also report similar behaviour in 
(2+l)-dimensional QED. Authors of this work employ 
the Schwinger proper time method, neglecting the field 
dependent phase of the fermion propagator. Therefore, 
a direct comparison with our findings is not straightfor- 
ward. Contrary to the widely studied case when the field 
is strong and the LLL dominates, all the Landau lev- 
els should be taken into account in the weak field limit. 
This feature makes the problem a difficult one and hence 
has been discussed less frequently in literature. Here we 
have shown that under plausible assumptions about the 
behavior of the mass function, the sum over Landau lev- 
els can be performed. The relaxation of some of the as- 
sumptions made is a natural generalization of this work, 
along with the inclusion of a thermal bath and the study 
of the gauge dependence of the results in the context 
of the Ward identities liSl and the Landau-Khalatnikov- 
Fradkin transformations . 
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